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What is claimed is: 

1. A computer implemented method for calculating a normalization factor 
comprising: 

providing a first intensity value (f } ) of a probe in a first probe array and a 
second intensity value (/ 2j ) of said probe in a second probe array; 

obtaining the geometric mean (x) of said / 7) and said / 2j ; 

calculating said normalization factor according to: 

f(x) - e hix) , wherein said h(x) is derived from referential intensities from 
said first and second probe arrays. 

2. The method of Claim 1 wherein said h(x) is derived by relating geometric means 
(xi') of first referential intensities (Rlf 0 ) in the first probe array and second 
referential intensities (RIi (2) ) in the second probe array to: 



, *A C2) , 

v 1 J 



3. The method of Claim 2 wherein said relating comprising: 

sorting (x it yd pairs according to Jt, into a plurality (m number) of bins with no 
overlapping; 

computing medians (x k ) of Xi's and medians ( y k ) of y t 's for each bin; and 
20 interpolating said medians ( x k , y k ). 
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4. The method of Claim 3 wherein said bins are of approximately equal size. 

5. The method of Claim 4 wherein said h(x) is: 

if x < x\ 

h{x) = 

if x > x 



jtj + 1 - x 



wy x + (1- w)y 1+p if xe (x u x\ +l),w =- ,/ = l,...,m-l, 

X j + 1 - X j 



Za 6. The method of Claim 5 wherein said m is 3. 

Hi 7. A computer implemented method for comparing the expression of a gene in a first 

Hi 

IS 10 sample with a second sample comprising: 

w providing a first plurality of intensity values ( ), each of which reflects 

the expression of said gene in said first sample, wherein said intensity values are 
obtained from a first nucleic acid probe array; 

providing a second plurality of intensity values (/ t (2) ), each of which 
15 reflects the expression of said gene in said second sample, wherein said intensity 

values are obtained from a second nucleic acid probe array; 

calculating a p-value using one-sided Wilcoxon's signed rank test, wherein 
the p-value is for a null hypothesis that median(f(x) l\ 2) - )=0 and an alternative 
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hypothesis that median((f(x) - l\ 2) )>0, wherein smdf(x) is a normalization 
factor; and 

indicating whether said transcript is present based upon said p-value. 

5 8. The method of Claim 7 further comprising a step of calculating normalization 
factor, said step comprising: 

obtaining the geometric mean (x) of said and said 7 ( (2) ; 

calculating said normalization factor according to: 

f(x) = e h{x) , wherein said h(x) is derived from referential intensities from 

10 said first and second probe arrays. 

9. The method of Claim 8 wherein said h(x) is derived by relating geometric means 
(jc/0 of first referential intensities (Rl! I} ) in said first probe array and said second 
referential intensities (RIi (2) ) in said second probe array to: 



15 y f =log 



RI 



(2) 

' J 



10. The method of Claim 9 wherein said relating comprising: 

sorting (jc,-, yd pairs according to x z into a plurality (m number) of bins with no 
overlapping; 

20 computing medians (x k )of jc,- 's and medians ( y k ) of y t 's for each bin; and 
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interpolating said medians (x k ,y k ). 



11. The method of Claim 10 wherein said bins are of approximately equal size. 

5 12. The method of Claim 1 1 wherein said h(x) is: 

if x < xi 

y L' - _ *< + 1 " x 
h(x) = i wy t + (1 - w)y i+1 , if xe (x i9 Xi + 1], w = -= — , / = l,...,m - 1, 

— ^ x i + i + ~ x i 

ifx>X m . 

13. The method of Claim 12 wherein said m is 3. 

14. A system for calculating a normalization factor comprising: 
10 a processor; and 

a memory coupled with the processor, the memory storing a plurality of 
machine instructions that cause the processor to perform a plurality of logical 
steps when implemented by the processor, the logical steps comprising: 

providing a first intensity value of a probe in a first probe array and a 
15 second intensity value (f 2) ) of said probe in a second probe array; 

obtaining the geometric mean (x) of said t }) and said f 2) ; 

calculating said normalization factor according to: 

f(x) - e h(x) 7 wherein said h(x) is derived from referential intensities from 
said first and second probe arrays. 
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15. The system of Claim 14 wherein said h(x) is derived by relating geometric means 
00 of first referential intensities (Rlf 1 *) in the first probe array and second 
referential intensities {Rlj 2) ) in the second probe array to: 



( Rl^ 



pr(2) 

v* 7 < J 



16. The system of Claim 15 wherein said relating comprising: 

sorting (x iy yi) pairs according to jc< into a plurality (m number) of bins with no 
overlapping; 

10 computing medians ( x k ) of x t 's and medians ( y k ) of 's for each bin; and 

interpolating said medians (x k , y k ). 



17. The system of Claim 16 wherein said bins are of approximately equal size. 



15 18. The system of Claim 17 wherein said h(x) is: 



h(x) = 



_ if x < x\ 

wy { + (l-w)y 1+1 , if xe (x u xi +l),w = 

if X ^ Xm « 



X j + 1 - JC 



X j + 1 - JC j 



,/ = l,...,m-l, 



19. The system of Claim 18 wherein said m is 3. 
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A system for comparing the expression of a gene in a first sample with a second 
sample comprising: 

a processor; and 

a memory coupled with the processor, the memory storing a plurality of 
machine instructions that cause the processor to perform a plurality of logical 
steps when implemented by the processor, the logical steps comprising: 

providing a first plurality of intensity values (7, (1) ), each of which reflects 
the expression of said gene in said first sample, wherein said intensity values are 
obtained from a first nucleic acid probe array; 

providing a second plurality of intensity values ( l\ 2) ), each of which 
reflects the expression of said gene in said second sample, wherein said intensity 
values are obtained from a second nucleic acid probe array; 

calculating a /?-value using one-sided Wilcoxon's signed rank test, wherein 
the p-value is for a null hypothesis that median(f(x) l\ 2) - )=0 and an alternative 
hypothesis that median((f(x) - 1\ 2) )>0, wherein said f(x) is a normalization 
factor; and 

indicating whether said transcript is present based upon said p-value. 

The system of Claim 20 further comprising a step of calculating normalization 
factor, said step comprising: 
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obtaining the geometric mean (x) of said / ( 0) and said // 2) ; 
calculating said normalization factor according to: 
f(x) = e h(x) , wherein said h(x) is derived from referential intensities from 
said first and second probe arrays. 

22. The system of Claim 21 wherein said h(x) is derived by relating geometric means 
(jc,-') of first referential intensities (RIi (1) ) in said first probe array and said second 
referential intensities (RIi (2) ) in said second probe array to: 

V 1 J 



23. The system of Claim 22 wherein said relating comprising: 

sorting (jc,-, yd pairs according to jc, into a plurality (m number) of bins with no 
overlapping; 

computing medians (x k )of xi's and medians ( y k ) of y t 's for each bin; and 
interpolating said medians ( x k , y k ). 



24. The system of Claim 23 wherein said bins are of approximately equal size. 



25. The system of Claim 24 wherein said h(x) is: 
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if x < x\ 



- X 



h(x) = < 



wy,+(l-w)y, 



if JCG (Xi,Xi + l],w = ^ 



,/ = !,. ..,m-l, 



y 



m ' 



if X > X m . 



The system of Claim 25 wherein said m is 3. 



A computer software product for calculating a normalization factor comprising: 



and said t 2) ; 

computer program code for calculating said normalization factor according 

to: 

fix) = e h{x) , wherein said h(x) is derived from referential intensities from 

said first and second probe arrays; and 

a computer readable medium for storing said codes. 

The computer software product of Claim 27 wherein said h(x) is derived by 
relating geometric means (*,-') of first referential intensities (RIi (1) ) in the first 
probe array and second referential intensities (Rlf 2) ) in the second probe array to: 



computer program code for providing a first intensity value (/ /; ) of a probe 
in a first probe array and a second intensity value (/ 2; ) of said probe in a second 



probe array; 



computer program code for obtaining the geometric mean (x) of said f 1} 
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y, =log 



\ 1 J 



111 



29. The computer software product of Claim 28 wherein said code for relating 
comprising: 

computer program code for sorting (*,-, yd pairs according to x, into a plurality (m 
number) of bins with no overlapping; 

computer program code for computing medians (x k ) of Xi's and medians (y k )of 



yi's for each bin; and 

computer program code for interpolating said medians (x k7 y k ). 



M= 10 



30. The computer software product of Claim 29 wherein said bins are of 
approximately equal size. 



31. The computer software product of Claim 30 wherein said h(x) is: 



15 



h(x) 



if x < xi 



wy x + (1 - w) y 1+1 , if xe (x x , xi + 1), w = 
ym ' ifx>x m . 



J = l,...,m-l, 



32. The computer software product of Claim 31 wherein said m is 3. 
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A computer software product for comparing the expression of a gene in a first 
sample with a second sample comprising: 

computer program code for providing a first plurality of intensity values 
( ), each of which reflects the expression of said gene in said first sample, 

wherein said intensity values are obtained from a first nucleic acid probe array; 

computer program code for providing a second plurality of intensity values 
( l\ 2) ), each of which reflects the expression of said gene in said second sample, 
wherein said intensity values are obtained from a second nucleic acid probe array; 

computer program code for calculating a p-value using one-sided 
Wilcoxon's signed rank test, wherein the p-value is for a null hypothesis that 
median(f(x) l\ 2) - 7 ( (1) )=0 and an alternative hypothesis that median((f(x) - 

l\ 2) )>0, wherein said/fx) is a normalization factor; 

computer program code for indicating whether said transcript is present 
based upon said /?-value; and 

a computer readable medium for storing said codes. 

The computer program code of Claim 33 further comprising computer program 
code for calculating normalization factor, said code comprising: 

code for obtaining the geometric mean (x) of said and said I- 2) ; 

code for calculating said normalization factor according to: 
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f(x) = e h(x) , wherein said h(x) is derived from referential intensities from 
said first and second probe arrays. 

35. The computer software product of Claim 34 wherein said h(x) is derived by 
5 relating geometric means (x,-') of first referential intensities (Rlf 0 ) in said first 

probe array and said second referential intensities (Rl! 2) ) in said second probe 
array to: 



v I 
V 1 J 



m 10 36. The computer software product of Claim 35 wherein said code for relating 
Uf comprising: 

u computer code for sorting (jc,-, yd pairs according to into a plurality (m number) 

of bins with no overlapping; 

computer code for computing medians ( x k ) of x ( 's and medians ( y k ) of for 

15 each bin; and 

computer code for interpolating said medians (x k ,y k ). 



37. The computer software product of Claim 36 wherein said bins are of 
approximately equal size. 

20 
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The computer software product of Claim 37 wherein said h(x) is: 



h(x) = 



if x < x\ 

i' _ _ _ x i + l ~ x 

wy { + (1 - w)y i+1 , if xe (xi,Xi + 1], w = ^ 

m if X ^ Xrn ■ 



Xi 4. i + 



,i = l,...,m-l, 



The computer software product of Claim 38 wherein said m is 3. 
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